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Abstract: At small x, the effects of finite transverse momenta of partons inside a hadron 
become increasingly important, especially in analyses of jets and heavy-quark production. 
These effects can be systematically accounted for in a formalism based on /cj_-factorization 
and unintegrated distribution functions. We present results for the unintegrated distribu- 
tion function, together with the corresponding integrated one, obtained within the frame- 
work of the Linked Dipole Chain model. Comparisons are made to results obtained within 
other approaches. 
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1. Introduction 

In the description of a given cross section in deeply inelastic lepton-hadron scattering 
(DIS), it is not enough to consider only the leading order perturbative terms. Although a s 
may be small, each power of a s may be accompanied by large logarithms due to the large 
phase space available for additional gluon radiation. It is, however, often possible to resum 
these emissions to all orders in the leading logarithmic approximation (LLA). 

The most familiar resummation strategy is based on DGLAP [|], ||, ||, Q evolution, 
which resums large logarithms of the virtual photon momentum transfer, Q 2 . Within 
this formalism, the cross-section for any given physics process is calculated using collinear 
factorization of the form 

a(x,Q 2 )=a (x,Q 2 )Y, [ — C a (z) f a {-,Q 2 ), (1.1) 

a •* 

convolution of coefficient functions C a and parton densities f a (x,Q ). 
DGLAP evolution describes most experimental results 1 from electron-proton and 
proton-proton colliders. By using input parton densities which are sufficiently singular 
when x — > 0, this formalism can also account for the strong rise of F2 for small x, as 
observed at HERA. However, there are problems with the description of non-inclusive ob- 
servables such as forward jet production in ep and heavy-quark production in ep and pp 
collisions. 

1 see e.g. pf for a recent review. 
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In the region of very small x (asymptotically large energies), effects of finite trans- 
verse momenta of the partons may become more and more important. The appropriate 
description in this region of phase space is BFKL evolution || |7) . The cross-sections are 
calculated in the /c_|_-factorization approach of the form 

a( x ,Q 2 ) = [ — d 2 k ± a(z,Q 2 ,kl)T(-,kl), (1.2) 
J z z 

convolution over the energy fraction z and the transverse momentum k± of the 
incoming parton of off-shell partonic cross-sections a, and A)j_-unintegrated parton densi- 
ties 2 , J-(x, This corresponds to a resummation of large logarithms of 1/x. The BFKL 
evolution equation actually predicts a strong power- like rise of F2 at small x. 

There exist a couple of models which take into account large logarithms of both Q 2 
and 1/x in DIS, reproducing both DGLAP and BFKL in the relevant limits. One such 
model, valid for both small and large x, has been developed by Ciafaloni, Catani, Fiorani 
and Marchesini, and is known as the CCFM model ||, [|. The resulting unintegrated 
distribution functions, J-(x,kj_,q), depend on two scales, the additional scale, q, being a 
variable related to the maximum angle allowed in the emission. 

The Linked Dipole Chain model (LDC) [JO], 11] is a reformulation and generalization of 
the CCFM model. Here, the unintegrated distribution functions are essentially single-scale 
dependent quantities, F(x, fc 2 ). In this article we present results for the integrated and 
unintegrated gluon distribution functions obtained within the LDC formalism and make 
comparisons with the CCFM model and with results from other formalisms. 

This article is organized as follows. We start in section ^ by giving a short introduction 
to the DGLAP and BFKL formalisms for deeply inelastic ep scattering, moving over to 
a description of the CCFM model. We end this section with a somewhat more elaborate 
description of the Linked Dipole Chain model for DIS, which is implemented in the Monte 
Carlo event generator LDCMC Jl2[ . This program can be used to describe both structure 
functions and exclusive properties of the hadronic final states. 

Since the gluon distribution functions are not experimental observables they are not 
uniquely defined, but depend on the formalism used. We discuss this problem, and some 
different approaches presented in the literature, in section ||| 

In section ||] we present our results for the unintegrated and integrated gluon distribu- 
tion functions, obtained in the LDC formalism as implemented in LDCMC. These results 
are also compared to those of other approaches, and we discuss how to make relevant 
comparisons between the different formalisms. 

We end this article with a summary in section 

2. Deep Inelastic Scattering 

Typically a deeply inelastic scattering event is represented by a fan diagram, as the one 
shown in figure [|. The (quasi-)real emitted gluons, constituting the initial-state radiation, 

2 In this paper we will use T for the unintegrated parton distributions in general, and Q for the uninte- 
grated gluon distribution, treated as densities in log i.e. Q{x) = xG(x). For the integrated ones we will 
use the standard notation / and g respectively. 
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are labeled qi, while the virtual propagators are referred to as fej. The figure represents an 
exclusive final state, with the final-state radiation explicitly marked as the dashed lines. 
The final-state emissions have to be defined so that they do not affect the cross-section 
and give negligible recoils to the emitting partons. The exact separation between initial- 
and final-state radiation depends, however, on the formalism used. This problem will be 
further discussed in section ||. 

We will here mostly discuss purely gluonic chains, which should give the dominating 
contributions at small x. We will use different approximations of the gluon splitting func- 
tion Pg S (z) = \ + ytr-j — 2 + z(l — z). Splitting a gluon means that we have two new gluons, 
carrying fractions z and 1 — z of the original gluon energy, which contribute to the gluon 
density, while the contribution from the original one must be subtracted. 

In analytical calculations the subtraction is achieved with the so-called plus prescrip- 
tion for the pole at z = 1 and the addition of a term proportional to 8(1 — z) to the 
splitting function. At asymptotically small x the leading contribution can be obtained by 
considering only the 1/z pole of the splitting function, thus effectively only adding one 
gluon, neglecting the recoil for the emitting mother gluon. In this way the problem with 
subtraction is avoided altogether. A third approach is to include both poles in the splitting 
function accounting for the two produced gluons. The subtraction is then handled by a 
Sudakov form factor, which multiplies each splitting and represents the probability that 
the gluon to be split has not already been split before. Here before has to be defined by 
the ordering imposed on the emissions (e.g. ^-ordering in DGLAP and angular order- 
ing in CCFM). In an approximation where the non-singular terms in P gg (i.e. the terms 
—2 + z(l — z)) are neglected, the Sudakov form factor, As, is given by 

/dq 2 olz 
©order, (2-1) 
qj_ 1- Z 

where a = 3a s /ir. This definition is used in the CCFM approach. Equivalently the 



Sudakov form factor can be obtained by considering energy- momentum conservation [13] 
and be written 

In A s = - J a-^zdzP gg (z)e oTder . (2.2) 

Here the non-singular terms in P gg can be included without problems, and this method is 
used in the LDCMC generator. 

To make the presentation more transparent, we will in the following discussion of the 
different approaches only use the asymptotically small x approximation, and will return to 
the Sudakov form factor in section ^. There we will also consider the convolution of the 
perturbative evolution with non-perturbative input parton densities, which will be ignored 
in the remainder of this section. 

2.1 DGLAP and BFKL 

In the DGLAP region, characterized by large Q 2 and limited 1/x, the dominating 
contributions come from -ordered chains which fulfill Q 2 kj_ n S> fej_ n _i ••■ an< ^ 
k + i > k+i+x- In the limit where also x is small, so that we can approximate the gluon 
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splitting function with P gg (z) ~ 1/z (in the double leading log approximation - DLLA) we 
can write the unintegrated gluon distribution function for a fixed coupling a s on the form 



„ n 



a- 



■0(xi-i -Xi)9(k 2 ± i -k\i_ x ) 
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In the case of a running coupling, a(Q 2 ) 
sion but with InQ 2 replaced by ln(lnQ 2 



and Xi = k +i /P +! tot- 
a /ln<5 2 



(2.3) 
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Qn+l 



we get the same exponential expres- 
) and a by «o- 

In the BFKL region of very small x and lim- 
ited Q 2 , chains that are not ordered in k± need 
to be accounted for, even though they are sup- 
pressed. The resulting unintegrated distribution 
function increases like a power at small x: 



1 



f(£u 



(2.4) 



with the function f(k±,x) describing a random 
walk in \n(k 2 ± /A 2 QCD ) 0, [H], Such a 

power-like behavior is in approximate agreement 
with HERA data, with A ~ 0.3. We note, however, 
that a corresponding increase is also obtained from 
NLO DGLAP evolution. 

Both the DGLAP and BFKL evolution were 
developed to describe inclusive quantities such as 
F2, but they can be interpreted as an explicit 
summation of initial-state bremsstrahlung (ISB) 
of quasi-real partons, and can thus be used to de- 
scribe exclusive multi-parton final states. To do 
this we must also include the final-state brems- 
strahlung (FSB) from the ISB partons within kine- 
matic regions allowed by the colour coherence constraint. This final-state radiation should 
also be emitted in such a way that it gives negligible recoils, and that it does not affect 
the total cross section. The separation between ISB and FSB depends upon the formalism 
used, and if more partons are treated as initial-state radiation we get a larger number of 
contributing chains, which is compensated by smaller weights for each chain, and with 
correspondingly reduced kinematic regions for final-state emissions. 

2.2 CCFM 

The particular calculation scheme adopted by Ciafaloni, Catani, Fiorani and Marchesini 
has resulted in the well-known CCFM model [||, Sj. This model, that has been developed 



proton 



Figure 1: A fan diagram for a DIS 
event. The quasi-real partons from the 
initial-state radiation are denoted qi, and 
the virtual propagators hi. The dashed 
lines denote final-state radiation. 
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assuming purely gluonic chains, provides a description not only of the structure function 
evolution in DIS, accurate at the leading-log level, but also of final-state partons. Colour 
coherence implies that the initial-state emissions are ordered in angle (or equivalently in 
rapidity). According to the definition of the separation between initial- and final-state 
radiation, they are also ordered in the positive (along the incoming proton) light-cone 
momentum q + . All other kinematically allowed emissions (symbolized by the q[ emission 
in figure |T|) are defined as final-state emissions. 

The CCFM model is based on the fej_-factorization formalism, with the unintegrated 
distribution function in the small- x limit given by: 



„ J z i 7r 9 I i 



. dzi d?q_[ 

S(x - ILz^efa - q i _ l z i - 1 )8{k'{ - k\ n )9{q - q n z n ). 



The notation is that of figure |], i.e. q±j and k±i are the transverse momenta of the real 
and virtual partons, respectively. The splitting parameter z is defined as z% = k+j/k+^-i, 
and the so-called rescaled transverse momentum q i is defined by q i = qj_i/(l — Zj). The 
interval for the z- variables is between and 1, which guarantees ordering in q+, and the 
angular ordering condition is satisfied by the constraint 

Qi>Qi-iZi-i, (2-6) 

explicitely written out in eq. ( |2.5| ). Moreover, we note the occurrence of the so-called non- 
eikonal form factor A ne (z, kj_, q) defined in ref. ||. The distribution function, Q, depends 
on two separate scales, the transverse momentum, k±, of the interacting gluon, and q, which 
determines an angle beyond which there is no (quasi-) real parton in the chain of initial- 
state radiation. In the rest frame of the incoming proton, this limiting angle corresponds 
to a rapidity given by (if counted negative in the direction of the probe) 

Mim = In {^Y) (2-7) 
In the original formulation there was also the so-called consistency constraint 

k 2 ±i > Zi q 2 ±i , (2.8) 
which was needed to guarantee that the virtuality k 2 is well approximated by — k\. This 



constraint has a non- leading effect, and has been disregarded in some analyses [17]. 

The CCFM evolution has been implemented in two hadron-level event generators, 
SMALLX [18, IE] and CASCADE pQ] , both maintained by Hannes Jung. These programs 



reproduce HERA data on F2 well for small x, where purely gluonic chains should give the 
dominating contribution. For larger x, we expect a large contribution from valence quarks; 
such chains are not easily accounted for in the CCFM formalism and are not included in 
the programs. 

Both programs are able to reproduce a wide range of final-state observables. We 
note, however, that there is a large sensitivity to non-leading corrections. In particular, 
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the description of forward-jet production at HERA turns out to be very sensitive to the 
non-singular terms in the gluon splitting function. In the original CCFM formulation these 
terms were left out, and without them the forward-jet rates are well described. If, however, 
the non-singular terms are included, which would be the most natural option, the jet rates 
come out approximately a factor two below the data |j. 

2.3 The Linked Dipole Chain Model 

The Linked Dipole Chain model (LDC) is based on the CCFM model, and agrees with 
CCFM to leading double log accuracy. Also LDC is formulated in terms of &_|_-factorization 
and unintegrated distribution functions. In LDC the ISB definition has been modified, 
resulting in a more simple description, with the unintegrated distribution functions being 
(essentially) dependent on only one scale and allowing for some sub-leading corrections to 
be introduced in a rather straight-forward manner. 

In LDC more gluons are treated as final-state radiation. The remaining initial-state 
gluons are ordered both in q + and g_ (which implies that they are also ordered in angle or 
rapidity y) with q±i satisfying 

q ±i > mm(k ±i , fcjL,i-i). (2.9) 

This redefinition of the ISB-FSB separation implies that one single chain in the LDC 
model corresponds to a set of CCFM chains. It turns out that when one considers the 
contributions from all chains of this set, with their corresponding non-eikonal form factors, 
they add up to one [10|. Thus, the non-eikonal form factors do not appear explicitly in 
LDC, resulting in a simpler form for the unintegrated distribution function 

/n ^2 
j J a — 1 - q +i )9{q-i - q- ti -i)5(x - Uzi) 8 (In k\ - lnk\ n ). 
z% nq Xi 

(2.10) 

The notation in the above and what will follow refers to that of figure |2[ Here, a typical 
DIS event is shown together with the corresponding phase space available in the 7-p rest 
frame, where the rapidity, y, and the transverse momentum, q±, of any final-state parton 
are limited by a triangular region in the (y, In g^)-plane. The proton direction is towards 
the right end of the triangle, and the photon direction is towards the left. The real emitted 
(ISB) gluons are represented by points in this diagram. The virtual propagators do not 
have well defined rapidities, and are represented by horizontal lines, the left and right ends 
of which have the coordinates (ln[k + i /k±i], In k^) and (ln[— In Ar^) respectively. 
The phase space available for FSB is given by the area below the horizontal lines (including 
the folds that stick out of the main triangle) . 

The ordering of the CCFM evolution in q + but not in q_ means that this formalism is 
not left-right symmetric. In contrast the LDC formulation is completely symmetric, which 
implies that the chain in figure [| can be thought of as evolved from either the photon or 
the proton end. (Thus, the LDC formalism automatically takes into account contributions 
from "resolved photons".) 
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In Q 2 bal/x 

Figure 2: The initial-state emissions qi in the (y,K = Zn(fc^))-plane. Final-state radiation is 
allowed in the region below the horizontal lines. The height of the horizontal lines determine hifc 2 ^. 
The light-cone momenta fc+i and k-i can be read off as described in the text. 



Returning to eq. ( |2.10| ), we note that it can equally well be expressed in terms of the vir- 



tual propagator momenta. Due to the condition in eq. ( |2 .9| ) we have q Xi ~ max(k^_ i , k 
and, suppressing the 6- and 5-functions, we obtain 

dzi 



dk 2 ±i 



(2.11) 



Zi max(fcj_ i ,fcj_ i _ 1 ) 

In particular we note that this implies that, for a "step up" or a "step down" in k±, the 
following weights result: 

d 2 a i i d 2 kj_i 



k 2 

K _Li 



k±i > kj_ j_i and 



d 2 q±i _ , d 2 k_ Li k] 
k 2 ' k 2 



k±i < kx 



t-i- 



(2.12) 
(2.13) 



i-1 



Thus, for a step down there occurs an additional suppression factor 

The relation between the integrated and the unintegrated distribution functions are 
symbolically written as xf(x,Q 2 ) ~ f ^^T(x, k 2 ^). The exact relationship is somewhat 
dependent on the evolution scheme used for J-(x, k 2 ^). For LDC we have 



xf(x,Q 2 ) 



dkj 
~kl 



■•F(z,fci) + 



dk 



•F (x ~ , & 



Q 2 



(2.14) 



assuming a single-scale dependent unintegrated distribution function. 



The first term of eq. ( 2.14 ) corresponds to chains whose struck parton is less virtual 



than the probe, k 2 , < Q 2 . The second term contains the suppressed contributions that come 
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from chains whose struck parton has k\ > Q 2 ; the suppression factor Q 2 jk\ is analogous 
to the one occurring in eq. ( |2.13| ), Note also that the x-argument in the unintegrated 
distribution function of the second term has been rescaled by the factor k^/Q 2 . 
Below we list some properties of the LDC model: 

• The natural scale in the running coupling, a s , is gj_-, which coincides with 
max (^L)^i,i-i)- 

• Non-leading effects such as those coming from quark- initiated chains, from the non- 
singular terms in the splitting functions and from energy-momentum conservation 
can be included in a straight-forward manner. It is also possible to correct the 
two emissions connected with the highest-virtuality link with the full 2 — ► 2 matrix 
element, thus improving the result. Finally, for the ISB emission closest to the photon 
end, the full off-shell 0{aa s ) matrix element is used. 

• The fact that fewer gluons are considered as initial-state radiation implies that typical 
z- values are smaller, thus resulting in smaller sub-leading corrections. 

• The formalism is greatly simplified by the fact that the non-eikonal form factors do 
not appear explicitly in the model. 

• The LDC formalism is fully left-right symmetric, in the sense that the result is in- 
dependent of whether the evolution starts at the probe (photon) or target (proton) 
end. Hence it can easily be generalized to the case when the photon is replaced by a 
hadron. As a result, the formalism can be applied to the study of jet production in 



hadronic collisions [21] 



Besides giving an inclusive description of the events, the result in eq. ( |2.11 ) can also 



be interpreted as the production probability for an exclusive final state. 

These qualities make LDC particularly suitable for implementation in an event genera- 
tor. One such LDC-based Monte Carlo event generator, LDCMC, has been developed by 



Kharraziha and Lonnblad [12]. 



The LDCMC has been shown to be able to reproduce available F2 data very well (see 
further section ||) , not only the small- x data from of HERA but also data at higher x from 
fixed target experiments. This success relies on the fact that quark- initiated chains are 
easily incorporated in the evolution. 

It should be mentioned that, just like the CCFM-based CASCADE, LDCMC has problems 
describing the HI forward-jet data. Again the data can be reproduced only when the non- 
singular terms of the gluon splitting function are omitted. The more physical case of 
including the full splitting function predicts a distribution which is approximately a factor 
2 below the data. 

In section || we present results for the integrated and unintegrated gluon distribution 
functions obtained within the LDC formalism using the LDCMC simulation program. Before 
comparing these results with those of other approaches, we discuss in the following section 
the relationship between the integrated and unintegrated distribution functions in different 
formalisms. 
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3. Discussion on Different Distribution Function Definitions 



In contrast to the structure function F2, the unintegrated parton distributions are not 
experimental observables, and depend on the specific scheme in which they are defined. 
Not only do they depend on the approximations used in the description of the evolution, 
it has also been argued that the unintegrated distributions are not gauge-invariant objects 
and there is a dependence on the gauge choice used for the off-shell matrix elements they 
are convoluted with 3 . 

For asymptotically large values of Q 2 or 1/x it is possible to define unintegrated par- 
ton distribution functions, ^(x,^), which depend on a single scale, k^_, specifying the 
transverse momentum (or virtuality) of the interacting parton. In the DGLAP region 
the fc^-ordering of the links implies that we can define J-(x, kj_) as the derivative of the 
integrated distribution functions, xf(x,Q 2 ), and thus 



dk 2 
k 



xf(x,Q 2 )= I ^T(x,k 2 ± ) (3.1) 



Also in the BFKL equation the distribution functions are determined by a single scale. 
However, since the chains are not ordered in k±, the integrated distribution function, 
f(x, Q 2 ), also obtains contributions from /cj_-values larger than Q 2 , and therefore the simple 
relation in eq. (EOI) is not satisfied. 



As discussed in section 2.2 the angular ordering in the CCFM model implies that 
the unintegrated distribution functions depend on two scales, k± and q, where q specifies 
the limiting angle in eq. (2.7). (Note that the strong ordering in k± in DGLAP, or in x 



in BFKL, also automatically implies an ordering in angle, which makes this second scale 
redundant.) In the CCFM formalism the gluon distribution function depends very strongly 
on the scale q, when q is in the neighbourhood of k±. The angular constraint is expressed 
by the last ^-function in eq. ( |2.5[) , which guarantees the relation 



q > q n z n , where q n = ^ n (3.2) 
This implies an upper limit for z n of the last emission given by 

z n < ^lim = (3-3) 

q + Q± 

From this result we can see that if the scale q is chosen in the neighborhood of k±, then a 
large fraction of the possible chains are cut away. To realize this we make two observations: 

(i) For q = k±_ we find z\i m ~ 0.5 when the last link is a step up in k± (in which case 
<1± ~ k±), and z\\ m < 0.5 when the last link is a step down (and q± > k±). 

(ii) In the CCFM model the 1/z pole in the splitting function is screened by the non- 
eikonal form factor, and the z-distribution obtained in the SMALLX and CASCADE MCs 
therefore does not peak at small z- values, but has a maximum around z = 0.5 B. 



3 See e.g. discussions in jB| and 
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This implies that for q = k± the constraint in eq. ( |3.3| ) will exclude a large fraction of the 
possible chains. Furthermore, the fact that for q = k±, zn m approximately coincides with 
the maximum in the z-distribution implies that increasing (decreasing) q in the neighbor- 
hood of k± includes (excludes) a significant set of chains. Consequently, for fixed k±, the 
structure functions depend strongly on q in this region. 

The relevant values for q in a hard sub-collision should, however, be significantly larger 
than k±. If the limiting angle is given by the final state parton in the hard collision, then 
it is easy to show that 

t = I (3.4) 
u 

where t, s, and u are the Mandelstam variables for the sub-collision. Thus if s is large 
compared to t we find q 2 ~ — t. For a very hard collision we find q 2 of the same order as s. 
We note that choosing q 2 = s corresponds to a limiting angle equal to 90° in the rest frame 
of the hard sub-collision. For a typical hard sub-collision we may thus have q substantially 
larger than k±, and we will return to this question in the following section. 

Many of the gluons which make up the initial-radiation chain in the CCFM model, 
are treated as final-state radiation in the LDC formalism. Therefore typical z-values are 
smaller in the LDC model, and most of the problem of angular ordering is postponed to the 
treatment of the final-state radiation. To leading order in lnl/x the result is determined 
by the 1/z pole, and the unintegrated distribution function in LDC depends on only a 
single scale, k\. As discussed in section |2[ sub-leading effects due to the 1/(1 — z) pole or 
the non-singular terms in the splitting function are included with Sudakov form factors, 
which do depend on the angular region allowed for radiation. Therefore also in LDC the 
unintegrated distribution functions depend on the scale q defined above, although as we 
will see below, the dependence is very much weaker for this model. 

Many schemes are presented in the literature to treat unintegrated parton distribu- 
tions. Besides with the CCFM formalism in the CASCADE and SMALLX MCs, which in 
the following will be referred to as JS (Jung and Salam) p3|, we compare in the following 
section our results also with the formalisms presented by Kwiecinski, Martin, and Stasto 
(KMS) Q and by Kimber, Martin, and Ryskin (KMR) (2J|. In KMS a term describing 
leading order DGLAP evolution is added to the BFKL equation. The parton distribution 
is described by a single scale, k±, and is assumed to satisfy the relation in eq. fl3,l|) . In 
KMR two-scale parton distributions are extracted from the same unified DGLAP-BFKL 



evolution equation, but as we discuss in more detail in section 4.2, the dependence on q 



for fixed k± is rather weak. Finally we will compare to a simple derivative, according to 



eq. (3.1), of the integrated gluon density from the GRV98 |2(| parameterization, referred 



to as dGRV in the following. 



4. Results 



In this section we discuss some results obtained from LDCMC. To illuminate the effect of 
the different contributions we study the following three different versions: 
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fit 
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b 9 
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&_L0 


fxg(x) 


X 2 /d.o.f. 
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1.86 
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1.78 


3 


0.57 


3 
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0.99 


0.37 


694/625 


gluonic 


2.71 





4 














1.80 


0.54 


193/86 


gluonic-2 


3.11 
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2.17 


0.39 


125/86 


leading 


2.34 





4 














1.95 


0.47 


126/86 



Table 1: The result of the fit of the parameters for the input parton densities. The standard 
version has been fitted to data from HI {27), ZEUS |2§, NMC |§] and E665 fo| in the region 
x < 0.3, Q 2 > 1.5 GeV 2 , while the gluonic and leading have been fitted to HI data only, in the 
region x > 0.013 and Q 2 > 3.5 GeV 2 . The last two columns give the resulting fraction of the proton 
momentum carried by the gluons and the x 2 over the number of fitted data points, respectively. 
Parameters in bold face have not been fitted. 



1. Standard: Including non- leading contributions from quarks and non-singular terms 
in the splitting functions. 

2. Gluonic: Including non-singular terms in the splitting functions, but no quark links 
in the evolution. 

3. Leading: Also purely gluonic chains, but with only the singular terms in the gluon 
splitting function. 

To get realistic results from these three versions we now also need to consider the convolu- 
tion of non-perturbative input parton densities. These are not a priori known, but need to 



be parameterized in some way. We will use the same parameterization as in 12] given by 

xfi(x,kl ) = A iX a *(l-x) b \ (4.1) 

where i = d v ,u v ,g and s for the d- valence, u- valence, gluon and sea-quark densities respec- 
tively (where the sea flavour densities are assumed to be fj=fu = 2/^). The parameters 
Ai,ai,bi and the perturbative cutoff, k±o, are then fitted to reproduce the measured data 
on i*2- There are some sum rules which fix the relationship between some of the param- 
eters. The A<i v and A Uv are fixed by flavour conservation and A s is fixed by momentum 
conservation. The fits to F% do not constrain the remaining parameters very strongly, so 
we have fixed the b parameters to 3 in the valence densities and to 4 in the sea and gluon 
densities. To check the sensitivity to the b parameter we have an additional fit for the glu- 
onic case with 6 = 7 called gluonic-2. In table [l] we present the result of the fits. Note that 
in the case of the gluonic and leading versions, only the gluon input density is considered. 

In figure || we show the resulting reproduction of F2 data. Clearly, all versions give 
a satisfactory fit to the data. The standard version gives an excellent description for all 
values of x. The gluonic and leading versions are naturally unable to fit the large- x data, 
as the inclusion of quark- initiated chains is essential for a description of large x- values. 
These versions are therefore only fitted to the small- x HI data, and these fits are shown 
separately in figure ^b. We note that gluonic-2 gives a better fit than gluonic, but we will 
anyway in the following concentrate on the latter, to make the comparison with the JS 
results (which also uses b g = 4) more informative. 



- 11 - 



(a) 




10 4 0.001 0.01 0.1 1 




x 



Figure 3: The description of F2 data as a function of x for the different fits presented in this paper. 
To separate the results for different Q 2 , what is shown is F2 + i for Q 2 . In (a) both small and large 
x data from HI, ZEUS @, NMC || and E665 §(| arc included. The full line is gluonic, dashed is 
leading and the dotted line is standard. In (b) the versions with only gluonic chains are compared 
to small- x HI data |27|]. As in (a), full line is the gluonic case, dashed is leading, whereas dotted is 
gluonic-2 (with b g =7). 

It should be noted that the version of LDCMC used here has not been released yet, and 
differs somewhat from the one described in Jl^] . The main difference is the handling of the 
Sudakov form factors, which were not quite correct in the original version. Also the full 
off-shell 7*5* — ► qq matrix element is now included 4 . These, and other minor changes, do 
not give a big effect on the results. 

In the following subsections we discuss the corresponding results for the integrated and 

4 Note, however, that the matrix element in J3l| is used, rather than the one from [^2) used in SMALLX 
and Cascade 
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the unintegrated gluon distribution functions. As we will see, the results are very similar 
for the standard and the gluonic version. The differences in the cascade can here be com- 
pensated by the adjustment of the input distribution functions and the perturbative cutoff, 
k±Q. For the leading version, without the non-singular terms in the splitting function, we 
find, however, that the resulting differences are small but not negligible. 

4.1 Results for the Integrated Gluon Distribution Function 

We start by studying the results obtained for the integrated gluon distribution function. 
As mentioned in section |2.3] (cf. eq. ( |2.14| )), in the LDC model the relation between the 
integrated and the unintegrated gluon distribution functions is as follows: 

xg(x, Q 2 ) = / -r±g(x, k\,Q) + ±g{ x ± k\ , Q)^- + xg (x, Q 2 ) x A s . (4.2) 

Thus, the integrated distribution function receives contributions from three different terms. 
The first term corresponds to struck gluons of transverse momenta below the virtuality of 
the probe, k 2 ^ < Q 2 (the full line in the triangular phase space of figure || is an example 
of such a chain), while the second term originates from chains whose struck gluons have 
k 2 ^ > Q 2 (e.g. the long-dashed line in figure^). Finally, the third term is the contribution 
from the input distribution function, corresponding to the case when no evolution has taken 
place. 



As discussed in section 2J3, the distribution function Q{x,k 2 _,Q) in the LDC model 
depends only very weakly upon the scale Q. This dependence is due to the Sudakov form 
factor corresponding mainly to the 1/(1 — z) term in the splitting function. Although this 
dependence is weak and can essentially be neglected in the first two terms in eq. ([L2]), 
it has a larger effect on the last term, where it is explicitely included as a multiplicative 
factor. This term dominates at large x, and the effect of the form factor is a suppression 
for larger values of Q 2 . 

We first study the relative importance of the different terms contributing to the LDC 
integrated gluon distribution function. In figure]^ the integrated gluon distribution function 
is shown as a function of x, for fixed Q 2 = 16 GeV 2 . As can be seen, at small x- values the 
first term of eq. ( [4.2j ) dominates; as x decreases the second term becomes noticeable. At 
large x-values, the behavior is governed by the input distribution function. 

Next we compare the LDC results with those obtained by other analyses. In figure |6| we 
show the LDC integrated gluon distribution function as a function of x for Q 2 = 16 GeV 2 
and Q 2 = 100 GeV 2 . We show our result for the gluonic together with the gluonic-2 case. 
As discussed above, they both give a good fit to F2. We see that the LDCMC results lie 
significantly below the JS curve for large x, but above JS for smaller x-values. Also shown 
in figure § are the corresponding results for CTEQ5M1 @ and MRST20011 @. An 
important point is that, while the LDC and JS results have been fitted to F% data only, the 
CTEQ and MRST curves have been fitted to more data. We see that these latter curves 
are more or less in agreement at large x, where there is more data available, while they 
separate more for smaller x-values. Clearly the LDC result agrees well with these curves 
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Figure 5: Relative importance of the different terms contributing to the LDC integrated gluon 
distribution function, for fixed Q 2 = 16GeV 2 . The dotted curve corresponds to keeping only the 
first term in eq. (4.2); keeping the first and second terms we obtain the dashed curve. The total 
result is represented by the full curve. 



for small x, where it almost coincides with CTEQ5M1. For larger x the LDC curve lies 
above CTEQ and MRST for gluonic, but agrees well with them for gluonic-2. 

This comparison leads to increased confidence in the physical relevance of the LDC 
model in general and in the LDC unintegrated gluon distribution function in particular. 
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Figure 6: The LDC integrated gluon distribution function (full curve is gluonic, dotted curve is 
gluonic-2), compared to the corresponding results of JS (dash-dotted curve), CTEQ (short-dashed 
curve) and MRST (long-dashed curve), for (a) Q 2 = 16 GeV 2 and (b) Q 2 = 100 GeV 2 . 
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Figure 7: The LDC unintcgrated gluon distribution function as a function of (a) x and (b) 
k\ . Study of sub- leading effects caused by inclusion of quarks and non-singular terms in the gluon 
splitting function. Full line is gluonic. dotted line is leading and dashed line is standard. (In these 
figures q is put equal to k±.) 



We shall study this in the next subsection. 

4.2 Results for the Unintegrated Gluon Distribution Function 

We now turn our attention to the topic of unintegrated gluon distribution functions, and 
we first want to study the effects from sub-leading corrections caused by the inclusion of 
quarks and by the inclusion of the non-singular terms in the gluon splitting function, in 
order to verify the statements made earlier in this section. 

In figure |7| the unintegrated LDC gluon distribution function is shown both as a func- 
tion of x and k± . Comparing the result for the purely gluonic case to that when we allow for 
quarks as well, we see that the differences are very small. (Note that the input distribution 
functions have been refitted.) Even though the effect of omitting the non-singular terms 
in the gluon splitting function is rather small, it is nevertheless not completely negligible: 
as can be seen in figure ^a, at small x there is a discrepancy of about a factor of two, a 
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Figure 8: The LDC gluonic unintcgratcd gluon distribution function (full curve), compared to 
the corresponding results of JS (long-dashed curve) and KMR (dotted curve), as functions of q/k± 
for (a) x = 0.01 and k± = 3GeV and (b) x = 0.001 and k±_ = 7GeV. 
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Figure 9: The LDC gluonic unintcgratcd gluon distribution function (full curve), compared to 
the corresponding results of JS (long-dashed curve) and KMR (dotted curve), as functions of k\ 
for (a) x — 0.001 and (b) x = 0.01. (In these figures q is put equal to 10GeT^.) 



result that is related to the problem with the description of HI forward-jet data that was 
mentioned in Section [2|. 

Our conclusion is that the non-leading effects can be largely compensated by slight 
modifications of the input distribution functions; we note in particular that the results 
with and without quark links are almost identical. 

We now want to study the dependence upon the two different scales discussed in section 
P| We show in figure || the dependence on q for fixed k± for the LDC, CCFM, and KMR 
formalisms. As discussed in section [3|, the g-dependence is rather weak in the LDC model, 
but very strong in the CCFM approach. We also see that the result in the KMR formalism 
is quite insensitive to variations in q. This is also illustrated in figure |||, which shows the 
results as functions of k\ for fixed q = 10 GeV. We see that the JS results start to fall 
dramatically below the other two, as k± approaches q. 

From figure || we note, however, that the CCFM result saturates for q > 2k±. In a 
hard- interaction event the scales |t| and s are normally larger than k\, and often charac- 
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Figure 10: The LDC gluonic unintegrated gluon distribution function (full curve), compared 
to the corresponding results of JS (long-dashed curve), KMR (dotted curve), KMS (short-dashed 
curve) and dGRV (dash-dotted curve) as functions of k\ for (a) x = 0.001 and (b) x — 0.01, and as 
functions of x for (c) k\ = 10 GeV 2 and (d) kj_ = SOGeV 2 . Results for LDC, JS and KMR, with 
q = 2k±, arc shown together with the 1-scaled distribution functions of KMS and dGRV. 



teristically by a factor of this order. For this reason we want to argue that when comparing 
the different formalisms, it is more relevant to study the CCFM distributions for q equal 
to the saturation value, rather than e.g. for q = k±. Indeed, for this larger q- value we see 
in figure | that there is a rather good agreement between the three models. 

This feature is further illustrated in figure 1C, which shows the distribution functions 
for q = 2k j_, as a functions of k 2 ^ for fixed x and as a functions x for fixed k±, and we see 
indeed a reasonable agreement between the LDC, JS, and KMR results. In these figures we 
also show the single scale KMS and dGRV results. Although these earlier parameterizations 
are somewhat lower for larger ^-values, we note a fair overall agreement between all five 
models. 



5. Summary 

Unintegrated parton distribution functions are not uniquely defined. They are not experi- 
mental observables, and their definition and properties depend critically on the formalism 
used. Thus when calculating e.g. cross-sections for production of jets or heavy quarks it 
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is necessary to use a consistent formalism for off-shell parton cross-sections and parton 
distribution functions. 

In this paper we present results for integrated and unintegrated gluon distribution 
functions according to the definitions in the Linked Dipole Chain model, obtained from the 
LDCMC program. We compare them with those obtained in other formalisms, in particular 
with results from the CCFM model obtained from the SMALLX and CASCADE MCs, and 
we demonstrate how to make a relevant comparison between the models. Indeed we find 
in this way a reasonable agreement between distributions obtained in different formalisms. 

Adjusting the input distribution functions for k± = k±Q, it is possible to find a good fit 
to the structure function F2 from LDCMC. The corresponding integrated gluon distribution 
function agrees well with fits to more complete data sets obtained by the CTEQ or MRST 
collaborations. The result is rather insensitive to non-singular contributions in the gluon 
splitting function or from from quark links in the chain. The contributions from quarks can 
be almost fully compensated by adjusting the input distribution functions. Omitting the 
non-singular terms in the gluon splitting function has a somewhat larger effect, reducing 
the gluon distribution for small x by roughly a factor 2. 

In the CCFM formalism the unintegrated parton distribution functions depend sen- 
sitively on two different scales: k±, which specifies the transverse momentum and the 
virtuality of the interacting parton, and q, which determines an angle beyond which there 
is no (quasi-)real parton in the chain of initial-state radiation. Many of the gluons which 
make up the initial-radiation chain in the CCFM model are treated as final-state radia- 
tion in the LDC formalism. Therefore typical z-values are smaller in LDC, and most of 
the problem of angular ordering is postponed to the treatment of the final-state radiation. 
This implies that in the LDC formalism the gluon distribution function is quite insensitive 
to the second scale q, and to leading log 1/x it depends only on a single scale k±. 

In a typical hard sub-collision the relevant additional scale q should be given by q 2 ~ \t\ 
or s. We observe that in the CCFM formalism, for a fixed k±, the gluon distribution 
function saturates for q > 2k±. This corresponds to typical values of t or s, and we therefore 
suggest that for a relevant comparison between the different formalisms we should choose 
q in this saturation region. Here we find indeed a reasonable agreement, not only between 
LDC and CCFM, but also with e.g. the two-scale formalism presented by Kimber, Martin 
and Ryskin, and the one-scale formalism by Kwiecinski, Martin and Sutton. Actually we 
see that it also agrees rather well with the derivative of the integrated distribution function 
GRV98, obtained in the pure DGLAP formalism. 
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